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HIGHER-ORDER VARIATIONAL CALCULUS ON LIE 

ALGEBROIDS 

EDUARDO MARTINEZ 


Abstract. The equations for the critical points of the action functional de¬ 
fined by a Lagrangian depending on higher-order derivatives of admissible 
curves on a Lie algebroid are found. The relation with Euler-Poincare and 
Lagrange Poincare type equations is studied. Reduction and reconstruction 
results for such systems are established. 


1. Introduction 

In this paper we study optimization problems defined by a cost functional which 
depends on higher-order derivatives of admissible curves on a Lie algebroid. Exam¬ 
ples of this type of problems are the optimal control of dynamical systems where 
the system to be controlled is a mechanical system, and hence depends on acceler¬ 
ations nmnuiH], trajectory planning problems in control theory [26] , key-framed 
animations in computer graphics m, and in general, problems of interpolation 
and approximation of curves on Riemannian manifolds [3 [27]. In many of these 
examples the presence of symmetries is used to reduce the difficulty of the problem. 

The advantage of the Lie algebroid approach is its inclusive nature, under the 
same formalism we can describe many systems which are apparently different m 
uniiniiii] and hence it allows a unified description even in the case of a reduced 
system when symmetries are present in the problem. An alternative approach con¬ 
sists in a case by case study using Euler-Poincare and Lagrange-Poincare reduction 
techniques as in [mill]. 

Previous work on the variational description of first-order Lagrangian systems 
defined on Lie algebroids provides a convenient departure point for the generaliza¬ 
tion presented here. In |36j it was shown that the Euler-Lagrange equations for a 
Lagrangian system defined on a Lie algebroid are the critical points for the action 
functional defined on an adequate Banach manifold of admissible curves satisfying 
boundary conditions. Such a manifold structure is a foliated one, the leaves being 
the E-homotopy classes of admissible curves. 

It is frequently argued that the variational principle for reduced systems, and in 
general for systems on Lie algebroids, is a constrained variational principle, in the 
sense that some additional constraints are imposed to the admissible variations. 
The point of view of [36] and the present paper is different. The set of admissible 
curves where the action is defined is endowed with a reasonable Banach manifold 
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structure. The tangent space to such manifold of admissible curves is precisely the 
whole set of infinitesimal admissible variations. Therefore no additional constraint 
to the infinitesimal variations is imposed, or in other words, they are as constrained 
as in the standard case when formulated directly in the tangent bundle instead of 
on the base manifold. In such standard case the relevant topology on the set 
C^{I,TM) is the one induced by the manifold structure of C^(]R,M), which is 
a foliated structure where each leaf (connected component) is a the set of curves 
of the form 7 with 7 in a given homotopy class. A similar construction can be 
performed in the case of admissible curves on a Lie algebroid E by using the notion 
of A-homotopy [131ES] . 

Whether the reader agrees with the above argument or not, it should be clear 
that the variational principle stated here is a bona fide standard variational prin¬ 
ciple, as it is needed in optimization problems: the solutions of the higher-order 
Euler-Lagrange equations are the critical points of the action functional which is a 
smooth function on a Banach manifold. In this sense one should notice that some 
generalized variational principles that appeared in the literature [20l El] do not 
satisfy this property. 


Description of the results and organization of the paper. For a Lie algebroid E 
we consider the set E^ of (fc — l)-jets of admissible curves on E. Given a func¬ 
tion L€C°°{E^), which will be called the Lagrangian, we want to find the max¬ 
ima/minima/critical points of a cost/action functional S dehned by 

rti 

S{a)= L{a^{t))dt 
Jto 

among the set of admissible curves which are E-homotopic to a given admissible 
curve oo and satisfy some boundary conditions. Finite variations as(f) = a(s, t) are 
given by E-homotopies and infinitesimal variations are just the s-derivative of Oj. 
E-homotopies are special morphisms of Lie algebroids = a{s,i)dt -I- fi{s,t)ds 
which satisfy fl{s,tQ) = /3(s,ti) = 0. This condition corresponds to the fixed 
endpoint condition in the standard case. The corresponding infinitesimal variation 
depends only on the values of cr(t) = /3(0,t) and are denoted Since each 

E-homotopy class is a Banach manifold [131ES] ; we can properly talk about critical 
points of the function S. 

The differential equations satisfied by the critical points are called the higher- 
order Euler-Lagrange equations, and generalize the first-order Euler-Lagrange equa¬ 
tions defined by Weinstein m (see also EH EH). It will be shown that, in par¬ 
ticular, this equations are the higher-order Euler-Poincare equations when the Lie 
algebroid is a Lie algebra [18] , the higher-order Lagrange-Poincare equations when 
the Lie algebroid is the Atiyah algebroid associated to a principal bundle EZ], or 
the higher-order Euler-Poincare equations with advected parameters when the Lie 
algebroid is an action algebroid min], in addition to the standard higher-order 
Euler-Lagrange equations when the Lie algebroid is the tangent bundle to a mani¬ 
fold. 

One of the advantages of the formalism of Lie algebroids is that morphisms of 
Lie algebroids can serve to relate Lie algebroids of the different types mentioned 
above. When two Lagrangians are related by a morphism of Lie algebroids, the 
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corresponding variational problems are also related and this allows to easily de¬ 
duce correspondences between the critical points of the associated action function¬ 
als, which amounts to reduction theorems showing, among other results, that the 
standard higher-order Euler-Lagrange equations for a higher-order left-invariant 
Lagrangian (with/without parameters) on a Lie group reduce to the higher-order 
Euler-Poincare equations of the reduced Lagrangian (with/without advected pa¬ 
rameters) on the Lie algebra, or that the standard higher-order Euler-Lagrange 
equations for a higher-order invariant Lagrangian on a principal bundle reduce to 
the higher-order Lagrange-Poincare equations of the reduced Lagrangian on the 
Atiyah algebroid. 

The paper is organized as follows. In Section [2] we will introduce the necessary 
preliminary results about higher-order tangent bundles and Lie algebroids, and we 
will fix some notation. In Section [3] we will define the space of jets of admissible 
curves and we will study its basic properties. In Section|3]we will find the properties 
and the local expression of the variational vector fields defined by E-homotopies, 
and its relation to the complete lift of a section of E. In order to find an intrinsic 
expression of the Euler-Lagrange equations, we will define in Section [5] two dif¬ 
ferential operators, the variational operator and the Cartan operator. This will 
be done by introducing the vertical endomorphism. In Section [5] we will find the 
critical points of the action functional in terms of the variational operator and we 
will deduce its coordinate expression, as well as those of the Cartan form and the 
Legendre transformation. Also a version of Noether’s theorem follows easily from 
the variational character of the equations. The relevant typical examples are given 
in Section [T] Finally, in Section [ 8 ] we will study the transformation properties of 
critical points under morphisms of Lie algebroids which readily amount to reduction 
results among the different kind of equations for different Lie algebroids. 


2. PRELIMINARIES 

2.1. Higher-order tangent bundles. Let M be a manifold. For a curve 7 : R —>■ 
M, defined in some open interval containing the origin in R., we denote by [ 7 ]^ = /q 7 
the /c-jet of 7 at 0 , which is said to be the fcth-order velocity of 7 or simply the 
fc-velocity of 7 . The set of fc-velocities of curves in M is a manifold T^M, known as 
the /cth-order tangent manifold to M. For fc = 1 we have T^M = TM, the tangent 
bundle to M. The projections : [ 7 ]^ i-A [ 7 ]^ define bundles 

T^M T'M M, for fc > ; > 0 . 


See [161114] for more information. 

A vector tangent to T^M can be described by a 1-parameter family of curves 
7 : M defined locally in a neighborhood of the origin in K^. Concretely, the 

family 7 s(t) = 7 (s,t) dehnes a curve in T^M by fixing s and taking /c-jets [ 7 ®]^. 


The vector ^[7s] 


s=0 


tangent to such curve at s = 0 is tangent to T^M at the 


point [ 70 ]^- With a different notation, that will be used in what follows in this 
paper, it is the vector [s >->• [t >->• 7 ( 5 , t)]^]^. 

For a curve 7 : R —>■ M we will denote by 7 ^^^ the curve 7 *^^^: R —>■ given 

by7('=)(t) = [sH>7(t-bs)]'’’. 
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There is a canonical injective immersion T{T^M), defined by 

= [s i-A [t i-A 7 (t + s)]^]^. Such an immersion allows to identify {k + 1)- 
velocities with the vectors tangent to T^M which are in the image of 

The canonical flip map on T^M will be denoted by Xk- T^TM —>■ TT^M. It 
can be easily defined in terms of 1-parameter families of curves by 

(2.1) xk{[t ^ [s ^ 7(s,t)]T) = [s^[t^ 

Fixing local coordinates cc* in M, we have an induced system of local coordinates 
{xy, j=0,...,k,oi T'^M given by xyy) = ^{0). 

2.2. Lie algebroids. A Lie algebroid structure on a vector bundle t: E ^ M 
is given by a vector bundle map p: E ^ TM over the identity in M, called the 
anchor, together with a Lie algebra structure on the C''^(M)-module of sections of 
E such that the compatibility condition [a, frf\ = {p{a)f)r] + /[cr, 77 ] is satisfied for 
every / G C°°(M) and every tr, 77 G Sec(E). See [28] for more information on Lie 
algebroids. 

In what concerns to Variational Calculus and Mechanics, it is convenient to think 
of a Lie algebroid as a generalization of the tangent bundle of M. One regards an 
element a of if as a generalized velocity, and the actual velocity v is obtained when 
applying the anchor to a, i.e., v = p{a). A curve a: [toj^i] —> if is said to be 
admissible, or an if-path, if x{t) = p{a{t)), where x{t) = T{a{t)) is the base curve. 

The Lie algebroid structure is equivalent to the existence of a degree 1 derivation, 
d: Sec(A^if*) —>■ Sec(A^''‘^if*), which is a cohomology operator = 0, and is 
known as the exterior differential on E. A morphism of Lie algebroids is a vector 
bundle map $: if —>■ E' such that ^*od = do^*. We may also define the Lie 
derivative with respect to a section cr of if as the operator d^-. Sec(A^’if*) 
Sec(A^if*) given by do- = ia o d + doia- Along this paper, the symbol d stands for 
the exterior differential on a Lie algebroid while the non-slanted symbol d stands 
for the standard exterior differential on a manifold. 

A local coordinate system (a;*), i = 1,...,?7- = dim(M), in the base manifold 
M and a local basis {Ca}, a = l,...,m = rank(if), of sections of E determine 
a local coordinate system (a;*, 7 /“) on if. The anchor and the bracket are locally 
determined by the structure functions p), and on M given by 

d 

(2.2) p(ea) = Pa-^ and [ea,ep] = Clp e^. 

The exterior differential d on the Lie algebroid is locally determined by 

(2.3) da;* = pj,,e“ and de'’' = —^C2pe°‘ A e^, 
where {e“} is the dual basis of {ea}. 

The if-tangent to a fibration. [2911301 dS]. Let if be a Lie algebroid over a man¬ 
ifold M and tt: P ^ M he a fibration. The if-tangent to P is the Lie algebroid 
Tp : T^P P whose fibre at p is the vector space 

y^P = { (5, v) G E^(^p) X TpP I p{b) = Tp7r(u) } , 

the anchor is p(6, v) = v and the bracket is determined by the bracket of projectable 
sections. We will use the redundant notation (p, 6, v) to denote the element (6, v) G 
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j'Ep The projection onto the second factor {p,b,v) M- 6 is a morphism of Lie 
algebroids and will be denoted hy Ttt: T^P E. 

Given local coordinates (a;*, u^) on P and a local basis {ca} of sections of E, we 
can define a local basis {Xa, V^} of sections of T^P by 


Xa{p) = (p,ea{Tr{p)),Pa-^ ) and Va{p) = (p^O^vip ),)■ 


Locally, the Lie brackets of the elements of the basis are 


(2.4) [Xo,,Xp] = ClpX^, [X^,Vb\=Q and [14,Vs] = 0, 
and, therefore, the exterior differential is determined by 

dx^=plX°‘, du^ = V'^, 

(2-5) 1 

dX'r = AXl^, dV^ = 0, 

where {d:’“, is the dual basis to {Aji, Vyi}. 


First-order variational vector fields. Within the context of Variational Calculus on 
Lie algebroids, a variation of an admissible curve a: [tojG] -A E is associated to a 
morphism of Lie algebroids q;(s, t)dt + I3{s, t)ds : TR^ ^ E such that a{t) = q;(0, t) 
and /3(s, to) = 0, /3(s, ti) = 0. The variational vector held (0, t) is determined by 
a{t) = /3(0,t) and it is denoted SaCr(t). In local coordinates, it has the expression 


( 2 . 6 ) Ea{a){t) = pMt))a-{t) 


dx^ 


ait) 




a(t) 


where a and cr have local expression a(t) = ( 7 *(t), a“(t)) and cr(t) = cr“(t)ea( 7 (t)). 
See O [36] . 

The variational vector held can also be dehned in terms of the canonical involu¬ 
tion of the bundle T^E. See [121131] for the details. 


3. Jets of admissible curves 

Consider a Lie algebroicfl (r: M,p, [ , ]). A curve a: /cR E is said to 

be an admissible curve in E, or an A-path, if it satishes poa = 7 , where 7 = roa is 
the base curve. 

Definition 3.1. For fceN, we denote by E^ the set of ik — l)-jets of admissible 
curves on E\ 

E^ = { [a\^~^GT^~^E I a is an admissible curve in A } . 

These spaces were introduced by Colombo and Martin de Diego in [9]. See 
also [24] for a simple exposition. 

Remark. Notice the grading E^ = E, E^CTE, and in general E^gT^~^E. The 
notation is suggested by the intended use in higher-order mechanics, where the 
elements in E are already considered as quasi-velocities (i.e. 1-quasi-velocities are 
in E^ = E). For instance, in the standard case E = TM, we have E^ = TM, 
E'^ = T'^M, etc. In our notation, plain indices will indicate the space where the 
object is dehned, while indices between parenthesis will indicate jet prolongation 
or the number of derivatives. 


^Nearly all the material in this section remains valid for an anchored vector bundle (t: E —t 
M,p) 
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To gain some intuition, it is convenient to describe the situation locally. Taking 
local coordinates on E, an admissible curve a{t) = (7*(t),a“(t)) is deter¬ 

mined by the function and the initial value 7*(0), since the function 7*(t) can 
be determined as the solution of the initial value problem = p]^{x)a°‘[t) with 
initial condition a;(0) = 7(0). Therefore, the {k — l)-jet of a{t) corresponds just to 
the {k — l)-jet of the function a“(t) together with the initial value 7*(0), the 0-jet 
of 7*(t). The natural coordinates of [a]^~^€T^~^E are given by 

x\o) = 7*(0), 

(3.1) = r = l,...,fe-l, 

where 'kj, are smooth functions depending also smoothly on p), and its partial 
derivatives up to order r — 1, obtained by taking total derivatives of the admissibility 
condition = p^a’^. 

Conversely, given a point (xq, yf,..., 2 /^)gR" x R'""™ the admissible curve given 
by a“(t) = solution 7*(t) of the initial value problem i* = 

p’^{x)a°‘{t), a;*(0) = Xq, is an admissible curve whose {k — l)-jet has coordinates as 
given in with = y“ for r = 1,..., fc - 1. 

It follows that E^ is a smooth submanifold of T^~^E with dimension dimi?^ = 
n + km, and that we can take a local coordinate system (a;*,2/“) of the form given 
above, x* = y“ = yf^_^y 

We will denote by Tk,i: E^ ^ E^ the restriction of the natural jet bundle projec¬ 
tion T^_i i_i : T^~^E —>■ T^~^E. Then it is straightforward to prove the following 
result. 

Proposition 3.2. The set E^ is a submanifold ofT^E. The dimension of E^ is 
dim(i?^) = dim(M) -|- A:-rank(i?). For k > I > 0 we have that ■ E^ E^ is a 
smooth fibre bundle. For k = I + 1 it is an affne bundle. 

See Section [3 bellow for the construction of E^ for some concrete examples of 
Lie algebroids. 

For an admissible curve a: K —>■ i? we will denote by a^’: K — E^ the natural 
jet-prolongation of a to E^, given by 

a^(t) = [s i-A a{s + t)]^~^. 

Notice that with the above notations a^(t) = 

Remark. The manifold can also be defined as a subset of the i?-tangent to 
E^ by the following inductive procedure. Starting with E^ = E, and once we have 
constructed E^, we define as follows. Consider the submersion Tk^ : E^ M 
and the i?-tangent E^E^ to E^ with respect to such projection. Then we define 

E^+^ = { Z€T^e^ I T§,iz) = rTk,k-i{z )}. 

For instance, E'^ = {{a,a,V) G T^E}, is the set of admissible elements, denoted 
by Adm(i?) in [29j. This construction allows to consider E^~^^ as a submanifold 
E^+^ CT^E^. However, it is preferable to consider E’^ as a separate manifold, and 
to define an embedding into T^E^ as described in the next paragraph. 
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Canonical inclusion into the i?-tangent. There exists a canonical injective immersion 
map ik.i : defined by 

7fc,i([a]") = (H"-\ [a]°, [s^[t^ a{s + 0]'=-^) , 

which does not depend on the representative i?-path a{t) of the /c-jet 

In terms of the natural jet-prolongation of admissible curves, the canonical im¬ 
mersion is given by 

/ \ 

(t)j. 

The canonical immersion can be considered as a section of E^ along Tk+i,k- 
We will use the symbol T = ik,i to denote such a section, that is 

(3.2) r([a]'=)=(^a'=(0),a(0),^(0)). 


For k = 1 the section T was already introduced in 

The associated derivation dp maps functions defined on E'^ to functions defined 
on and will be called the total time derivative operator. The reason for that 

name is that for every function E&C°°{E^) and every admissible curve a we have 


More generally 


dTE{a^+\t)) = -E{a\t)). 




Associated to the coordinate system (a;by“) in E^ we have the local basis 
{Xa,Va} of sections of T^E^ E^, as in Section[2l 


d 

Xa{A) = ( 




m = Tkfi{A). 


To write more compact expressions we will use the notation = Xa. 
In local coordinates we have the local expressions 


(3.3) 
and 

(3.4) 


k k 

T = ^ ^ V^, 

i=i i=o 


, ^ ^ c.dE ^ ^ dE 


When if is a Lie algebroid, the operator dx acts also on p-forms on E^ producing 
a p-form on For the dual basis {A", V“} of the local basis {Xa, V))} we have 

. 35 . dTX‘^ = v?-c^y,x\ 

dTV“=V“+D r = l,...,fc 

which follows easily from the definition dx = doix -I- ix^d. 
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Jet prolongation of admissible maps. Consider a second Lie algebroid r': E' ^ M' 
with anchor p'. A vector bundle map $: E ^ E' is said to be admissible if it maps 
admissible curves into admissible curves. If M —> M' is the base map, then $ 
is admissible if and only if p'o^ = Tipop. 

An admissible map induces a map : E^ E'^ by jet-prolongation as follows: 
if a is an admissible curve then 3>oa is admissible and we set 

In other words, is the restriction of the {k — l)-jet extension : T^~'^E —>• 

Tk-iE' of $ to which takes values in E'^. It is clear from the definition that 
for /c > / > 0, which is also valid for fc = 0 if we set = tp. 

When we look at E^ as a submanifold of we have that is the re¬ 
striction to E^ of the map E^~^ —>■ E'^~^ given by {A,b,V) M- 

In other words, we have i which 

in terms of the operator T reads E'^^^oT = 


4. Variational vector fields and complete lifts 


As it was explained above, a vector tangent to the fc-tangent space to a manifold 
is determined by a 1-parameter family of curves. Accordingly, a vector tangent to 
E^ is determined by a 1-parameter family a{s,t) of admissible curves in E, by the 
same procedure: [s i-)- [t i-a is a vector tangent to E'^ at the point 

[t ^ a( 0 ,f)]'=-ieAV 

In the calculus of variations on Lie algebroids the families a{s,t) of admissible 
curves are given by morphisms of Lie algebroids (p: TR^ —>■ A of the form a{sE)dt + 
/3{s,t)ds. The variational vector field can be determined in terms of 

a{t) = (3{0,t) and its derivatives up to order k. It is a vector held along a^{t), 
where a{t) = a( 0 ,t), and will be denoted E^a{t). 

Indeed, let us hnd the local expression of the variational vector held associated to 
the morphism adt + pids. In local coordinates, the family a is ( 7 *(s, t), a^{s, t)) and 
the family f] is ( 7 *(s, t),/3^(s, t)). The fact that adt + fids is a morphism amounts 
to 


(4.1) 


dE 


= 


dE 

ds 


da^ 

~di 






where the local structure functions p^ and Cjf^ are evaluated at the point 7 ( 3 , t). 
The curve as^{t) in E^ dehned by the family a is given by 


x"=E{s,t), yi=a^{s,t), = —(s,t),..., ^-(s,t), 


and the coordinates of := -^as^{t) 


are 


s=0 


Taking into account the equations (14.11) we have 
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and 


V 




^ dt 
and hence we have 


= ^(0,0 + + 


V: 


dt^ 

It conclusion, we have 


^ ^- tK + 

rl+l —1 ^ ^ rliT—1 L ' I'-f J’ 


dt 


= 2,...,k. 


(4.2) 


"fe 


(T = P„CT 


,±_ 

dx^ 


^ f;La I fya „/3 71 ^ 


dy^' 


It follows that 5^ cr is a differential operator in a of order k (depends on [cr]*^*^)) 
and a differential operator in a of order k — 1 (depends on a^{t) = [/i H> a{t+h)]^~^). 

Remark. In the classical notation of the calculus of variations we have 


Sy^ =&°‘ + Cp^a^a'>, 5y^ =—5y'^_^, for r = 2 ,..., fc 


as it should be expected. 

The variational vector field S^tr can alternatively be defined directly in terms of 
(T and a, without reference to the morphism </>, by jet-prolongation as follows. 


Proposition 4.1. Let a he an admissible curve and let a be a section of E along 
TOO. Consider the associated first order variational vector field S^ct: R —^ TE, 
and its {k — C)-jet prolongation {'E.aCr)^^~^\t)GT^~^TE. Then, the vector field 
Xfc-io(SaCr)('="^) is tangent to E^ and coincides with 5^ct. 


Proof. It is enough to prove the proposition for t = 0; for t = we can just 
consider the curves d{t) = a(to + t) and aft) = afto + t). 

We take a{s,t) a family of admissible curves such that a(0,t) = a{t) and find 
a complementary fi{s,t) such that adt + jdds: TR^ —>■ i5 is a morphism of Lie 
algebroids with /3(0,t) = (j{t). On one hand, by construction, we have that [s i-)- 
\t 1 -^ a(s,t)]^“^]^ takes values in TE^. On the other hand, ^(0,t) = EaCrft), from 
where 



s=0 


[s^ [t^a(s,t)]^-V 


= X,_i([t^[s^a(s,t)]V-') 

= Xk-l{[t EaCr{t)f~'^) 


= {Xfc_io(S,a)('=-i)}(0), 


and the result follows. 


□ 


As an immediate consequence of the above proposition we have 
(4.3) TrkjoE'^a = Efa, for fc > / > 0. 

It follows that 5^(7 projects to the vector held p{a), i.e. TT^poE^a = p{cr), and 
hence it make sense the following dehnition. 

Definition 4.2. Let a be an admissible curve and cr be a section along 7 = too. 
The section a’f of T^E^ along given by 

is said to be the fcth-order complete lift of a with respect to a. 
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It is easy to see that every element in is of the form cr^’(O) for some 

section a. 


Complete lift of a section of E. The concept of variational vector field is related to 
the concept of complete lift of a section of E. 

Let 77 be a section of E and ^s) its flow [T31I2S]- The map $ 8 ^ : E^ E^ is a 
flow in E^ which defines a vector field X^gTE’^. This vector field is rfc^;-projectable 
over Xlj for k > I > 0, and Tfcp-projectable over p{r]). Therefore it allows to define 
a section G Sec {T^E'^) by 

rj\A) = iA,r,{rk,oiA)),X^iA)), AgEK 

The section rj^ will be called the /cth-order complete lift of the section rj. For k = 1 
the section coincides with the complete lift rj^ of the section rj as defined in [ 2 ^. 
From the definition it is clear that TTk,iori^ = rfoTk^u for k > />0, where rf = r] 
should be understood for 1 = 0 . 


Proposition 4.3. For any section r]G Sec{E) we have the following property 
(4.4) dToi^k =i^k+iodj'. 

Conversely, if Z is a section of E^~^^ which projects to T]GSec{E) and satisfies 
dxoi^k = izodx, then Z = . 

Proof. It is enough to prove the proposition over basic forms 9€Sec{E*), and 
over forms of the type d'^{dF) for F€C^{E) and r = 0,..., fc — I (we have omitted 
the pullbacks for simplicity), since they generate the set of sections of (T^E^)*. 
For 9g Sec{E*), in local coordinates if 0 = 6ae°‘ and 77 = ri°‘ea, we have 

dT9 = {dTeof)X^ + 0„(V“ - 

and 

77 I = 77 “ Ala + (dT??“ + Cpyr]'^)Va, 

from where 

{drO ,ri^) = {dTeJv° + ^'a(dr?7“) = dT{0,r]). 

We will prove that for any function FgC°°{E) and r = 0,... ,k we have 

dT{dlrdF,r]^+^) = (, 77 ’'+^ ), 


which can be equivalently stated in the form 


dj’d^r+idyF’ = d^r-+2d^^F’. 

Consider an arbitrary point A = and the flow of the section rj. 

Then the left hand side evaluated at A is equal to 


(dj- d^r+i dip E'ji^A) 


y{d^r+kdlpF)oay{0) 

^^[dlpFoy^oayit) 

dt ds 


[d^Fo($,oa)’’+^](t) 


_d _d 
dt ds 

d d r d*" ,,, , 
d if'*' ,, , 


s=0 


s=0 






t=0 


s=0 




s=0 
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and the right hand side evaluated at A is 


= -(4+iFo$-+2)(Al) 

' as s=o 

= ^(4+iFo$^+2oa''+2)(0) 


'’'+Vo($,oa)’’+ 2 )( 0 ) 

d 




ds dr+i 


{Fo^sOa){t) 


t=o 


s=0 


s=0 


s=0 


and both expressions are equal. 

Conversely, let Z be a section of such that Trk+ifioZ = -qaTk+ifi and 

izdx)^ = dxirjkX for every section A of {E^E^)*. In particular for A = dl^^S with 
r = 1,.. . ,fc we have iz{d^6) = dTirif^dlf^d = i^k+i{d^9), so that iz-rj^^+^dXpO = 0 , 
for r = 1,..., fc. Moreover, since Z projects to rj it follows that this last relation 
holds also for r = 0. Therefore Z — 77 ^+^ = 0. □ 

Notice that as a consequence of the above property we have that 

(4.5) dxodrjk = d^k+iodx, 
which follows from dr^r = doi^r i^rod and dodx = dx^d. 

Relation of complete lifts with variational vector fields. The relation between com¬ 
plete lifts of sections and variational vector fields is as follows. 

Let a be an admissible curve over 7 . Consider a section rj of E and define the 
section a along 7 by a{t) = 77 ( 7 ( 1 )). Consider on one hand the complete lift ctJ of 
cr with respect to a, and on the other the complete lift 77 ^ of the section 77 . Then 

(4.6) at = 

Indeed, it is shown in |36j that the section 77 and the admissible curve a define a 
morphism (f) = a{s, t)dt + /3(s, t)ds = ^s{a{t))dt + r]{ips{'j{t)))ds. At s = 0 we have 
Q;(0,t) = a(t), /3(0,t) = 77(7(t)) = a{t) and the vector field (t) = E^a is 

equal to A^(a^(t)). Therefore 

at{t) = {a\t),a{tlEta{t)) = Xt{a\t)) = {t)), 

form where at{t) = ri^{a^{t)) follows. 

Conversely, given a{t) along 7 (f) we can find a time-dependent section 77 such 
that = a{t). Using the construction for the time-dependent section as 

in [36] we also have a^{t) = rf{t,a^{t)). 

It follows that if the coordinate expression of 77 is 77 = 77 “ ea then the coordinate 
expression of the complete lift is 

k 

+ E dfr-^ldT^ + K- 


r—1 


( 4 . 7 ) 

Siinmilarly 

( 4 . 8 ) al{t) = a°‘(t)Xa{a^{t)) + 
where a = a°'{t)ea{'j{t)). 


Jr—1 


k 

r—1 


[a- +Ciya^]{tma\t)). 
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5. The vertical endomorphism and the variational operator 

Working with forms and their time derivatives, it is usefnl to have an operator 
which formally resembles integration, that is, it is a kind of inverse of dx- This 
is the role of the vertical endomorphism, which is an endomorphism of the vector 
bundle T^E^. We introduce the vertical endomorphism as an auxiliary tool to 
define two differential operators, the variational operator and the Cartan operator, 
which will be fundamental in the development of the calculus of variations. 

5.1. Vertical endomorphism. We will define the vertical endomorphism S by 
its action on sections of the dual bundle. We will make no notational distinction 
between S and its dual S*, both will be written S. 

The set of sections of E^)* is generated (with coefficients in C°°{E^)) by 
sections of the form {TTk,r+i)* d'j.O for sections 9GSec{E*) and r = 0,... ,k. In 
particular, if {e“} is a local basis of Sec{E*), then a local basis of Sec{{T^E'^)*) 
is {(T'rfc,r+i)*d^e“} for r = 0,..., A. Therefore, it is sufficient to give the action 
of S over sections of such a form. To simplify the notation, in what follows we will 
omit the pullbacks 

We define the vertical endomorphism as the (1,1) tensor held S: E^E^ T^E^ 
which satishes 

(5.1) S{dlj.6) = rd'E^e, 

where S{0) = 0 should be understood for r = 0. Notice also that for r = 1 we have 
Sidxe) = 6. 

We have to check the consistency of the above formula when taking linear com¬ 
binations of sections. Obviously there is no problem when taking sums, so that we 
have to check consistency when multiplying by functions on the base. If we take 
/6», with then 

j=o 

so that 

5(d5^(/0)) = E('^)(4“V)-5(4^) 

1=0 

1=0 

= E 

1=1 ^ 

= rdlf\.fe), 

consistently with the given dehnition. 

A section A of {T^E^)* is said to be Tfc^r-semibasic if it vanishes over elements 
in the kernel of Erk^r- In other words, A is Tfc^r-semibasic if there exists a map 
A: E^E'^ {E^E^y such that X(Z) = {X,ETk,r{Z)) for every Z&E^E^ In 
particular, A is Tfe^o-semibasic if there exists a map A: E^E^ —i E* such that 
r]{Z) = (A ,ETkfi{Z )) for every ZET^E'^. It is clear that the space of Tfe^j-semibasic 
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forms is generated (over C°°{E^y) by the family of sections of the form dJj.0 for 

r = 0 ,... ,L 

Proposition 5.1. A section A € Sec((T'®i?^)*) is Tk^i-semibasic if and only if 
5'*+! ( 77 ) = 0. The image of is the set of Tk^k-i-semibasic forms. In particular, 
we have that = 0 . 


Proof. For r = 0,... ,k, a simple calculation shows that 


SPidtj-e) 




if p > r, 
if p<r. 


It follows that SP{\) = 0 if and only A can be written as a linear combination (with 
coefficients in C°°{E^)) of forms of the type dtj<6 with r = 0,... ,p — 1. Therefore 
SP(X) = 0 if and only if A is rfc_p_i-semibasic. Moreover, from this expression we 
have that the image of SP are Tfe_fe_p-semibasic forms. A simple argument shows 
that every rfc_fc_p-semibasic form is in the image of S'^. □ 


Proposition 5.2. For any section AGSec((T^A^ ^)*) we have 
(5.2) S{dTX)-dT{SX) = X. 

Proof. It is sufficient to prove the proposition for a section A of the form A = dlj,9 
for r = 0,...,fc — 1. Using the definition of S we have 

SidxX) - driSX) = S{dljf^e) - dT{S{dlrO)) 

= (r + I)dy0 — dTird^O) = d^O = X, 

which proves the statement. □ 


Coordinate expression. Even though it is not needed in this paper, we will find 
the local expression of the vertical endomorphism. We take a local basis {ca} of 
sections of E and the dual basis {e“} of sections of E*. Then we have a local 
basis {Xa,Vf,} of sections of T^E^ and the dual basis {A’“,V“} of sections of 
{T^E^Y. We have to find the image by S of such sections. From the definition 
it follows that S{X°‘) = 0. To find the expression of 5'(Vf) we use equation (13.51) . 
dj-X^ = Vf — Cp^.^y^X'^, and thus 

S'(Vi“) = S{dTX°‘ + = S{dTX°‘) = dTS{X°‘) + X^ = A'“. 

Finally, in order to find S'(V“) for r>2 we recall equation ()3.5I1 . V“ = dTV“_i, and 
using Proposition 15.21 we have 

5(V“) = S(dTV^_Y = dT5(V“_i) + V“_i 


for r>2. It easily follows by induction that 

.5(V“) = (r - l)V“_i + 4”'5(Vr), r> 2 . 

Therefore we found 

k 

A = ^ [(r - l)v; 0 V “_1 + v; 0 d^A“] . 

r=2 


(5.3) 
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Remark. A more explicit local expression for S can be obtained as follows. The 
value of <S'(V“) can be written in the form 

5(V“) = rV“_i + [d^T^S{V^) - V“_i] r>2. 

The expression d^X°‘ — V“ is a (r — l)th total time derivative of a semibasic 1-form 

- V“ = 

which follows from dTX°‘ = V“ — . Therefore we can also write 

k 

(5.4) 5 = 0 ^ [rV; 0 V“_i - 0 . 

r^2 

In more compact way 

k k 

(5.5) 5 = ^ r V; 0 V“_i - ^ V; 0 dlf^C^y.X^). 

r—1 r—2 

5.2. Variational and Cartan operators. Having in mind the procedure of in¬ 
tegration by parts that will be needed in the calculus of variations, we can define 
two differential operators (see [4] for the standard case): the variational operator 
E, mapping sections of [T^E^)* into sections of given by 

(5.6) E(A) = A - dT{S{\)) + ^4(^^(A)) + ... + (-1)"^4(5'=(A)), 

and the Cartan operator § which maps sections of (T^E^)* into sections of (j~Ej^ 2 k-i 
given by 

(5.7) S(A) = 5(A) - |dr(52(A)) + ... + (-l)'=-il4-i(5'=(A)), 

for A G Sec{{T^E^)*). From the definition of E and § it is clear that 

(5.8) E(A) = A-dT(S(A)). 

Moreover, using the Proposition 15.21 a long but straightforward calculation shows 
that 5(E(A)) = 0, so that E(A) is r 2 fe,o-semibasic, and also 5^(S(A)) = 0, so that 
S(A) is r 2 fc-i,fe-i-semibasic. Moreover S(c?tA) = A. 

Remark. In a more systematic way, one can proceed as in [1], where a family 
of operators was introduced to study complete lifts of vector fields and other 
related properties. For s = 0,..., fc we define the operator D^. mapping sections of 
(T^E^Y into sections of Yf^E^^~^Y be means of 

k 

j=r 

for A G Sec((T^E^)*). Following the arguments in [i], it is easy to see that the 
image of are r 2 fe_r,fc-r-semibasic sections, and that we have the relations 

(5.9) D^-i = - droDr 

r! 

(5.10) 5oD^ = rD^+i. 

From the second relation (15.101) it follows that only Do and Di are relevant, while the 
others can be determined inductively by Dj.+i = ^5 o D,. for r>l. The variational 
operator is E = Dg and the Cartan operator is § = Di. Also from (15.101) for r = 0 
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we have S' o Dq = 0 so that E(A) is T2fc,o-semibasic. Applying ^ to (15.1011 for 
r = 1 we get 

S'=(Di(A)) = {k- l)!S(Dfe(A)) = is'=+i(A) = 0, 
so that §(A) is T2fe-i,fc-i-semibasic. 


6. VARIATIONAL CALCULUS 

Let J = and fix two points and Given a La- 

grangian function L G C°°{E^) we consider the action functional 

fti 

(6.1) S{a)= / L{a^{t))dt 

Jto 

restricted to admissible curves a in A such that a^“^(to) = Aq and = Ai. 

We look for critical points of such a functional. 

Of course, to speak about critical points, local maxima or local minima, has 
no meaning if we do not give explicitly a differential manifold structure to the set 
of such curves. In the case fc = 1 it was shown in [35] that the relevant Banach 
manifold structure is the foliated one P(J, E), where each connected component (a 
leaf) is an A-homotopy class of admissible curves. Therefore, in the higher-order 
case we will use the same structure and we impose the additional conditions coming 
from the boundary conditions. In an alternative but equivalent way, we can restrict 
S to an A-homotopy leaf (which is a Banach submanifold) and we use differential 
calculus to find the critical points. It follows that finite variations are those defined 
by A-homotopies with the given boundary conditions. 

We denote by mo, toiGM the base points mg = Tfe_i o(Ao) and mi = Tfe_i o(Ai). 
From the results in [35], the following result is easy to prove. 


Theorem 6.1. The set 


(6.2) V{J,E)^l = {a&V{J,E) \ a is and = Aq, = A^} 

is a Banach submanifold of E)'!fff. The tangent space to 'P{J,E)'^ at a point 

aGV{J,E)il IS 


(6.3) 


TaVUE)il = {E^,a 


a is and (to) = 0, (ti) = 0 | . 


Remark. If a(TP {J , E)'fff is an admissible curve, the connected component of 
a in 'P(J, A)™! is the equivalence class Ed a of admissible curves in E which are 
A-homotopic to a. Similarly, if aGP(J, A)^!), the connected component of a in 
V{J,E)^^^ is the equivalence class H'^cHa of admissible curves in E which are 
A-homotopic to a and have higher-order contact with a at the endpoints. 


It follows that the infinitesimal variations to be used are of the form with 
[o']^~^{ti) = 0, i = 0,1. Using the classical notation of the calculus of variations 

Sx^=p>^a‘^, 6y‘f=&‘^+C^y^a\ Sy^ =r = 2,...,k, 

with y^{ti) = 0 for r = 0,..., fc — I, i = 0,1. 


^In what follows the symbol S stands for the action functional. No confusion with the vertical 
endomorphism is possible since it will not be used explicitly. We will only need the operators E 
and §. 
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Given an admissible curve a we take a curve as on with ao = a, and 

the corresponding _B-homotopy a{s,t)dt + j3{s,t)ds where as{t) = a(s,t). Taking 
the derivative at s = 0, 



s=0 



d j-f 


'im 


s=0 


dt = f (dL(a^(t)),-^as^(t) ) dt. 

Jtn ds s=0 


Defining a{t) = P{0,t) we have that -^as^{t) = 5^CT(t) and hence 


s=0 


(6.4) (d5'(a) ,5acr) =-^S'(as) =/ (dL(a'=(t)), E;^cr(t)) dt. 

ds s=o 


to 


In the above expressions d stands for the standard exterior differential on a mani¬ 
fold. Using the exterior differential d of the algebroid and taking into account 

that a^{t) = (a^(t), cr(t), 5^tT(t)) we have 

(6.5) d5'(a)(SaCr) = [ {dL{a^{t)) dt. 

Jto 


Let us consider the sections E(dL) and Ol = S{dL), which are related by E(dL) = 
dL — drdL- We have 


{dL{a\t)),al{t)) = {ndL){a^\t)),al\t)) + ((dT0L)(a^'=(t)), af (t)) 

= (E(dL)(a2^t)),af(t)) + |(0L(a^'=-'(t)),af-i(t)). 

Taking into account that E(dL) is Tfe^o-semibasic we will consider the associated 
map SL: E^^ —>■ E*. Similarly, taking into account that 9l = S{dL) is T 2 k-i.k-i- 
semibasic we will consider the associated map El ■ —)• {T^E^~^)*. We then 

have 

( dL{a>^{t ), a': it)) = { 5L{a^\t )), a{t )) + | ( ELia^'^-^t )), (t)). 

Inserting this into the variation of the action we get 


dS'(a)(5oCT) = 


'to 


( 5L{a^\t )), a{t)) + it )), ) 


dt 


= / {5L{a^\t)),a{t))dt + {EL[a^^-\t)),a^,-\t)) 


After imposing the boundary conditions we arrive to 


( 6 . 6 ) 


dS{a){Ea<y) = (SLia^'^it) ,a{t)) dt. 


'to 


Since cr is arbitrary, from the fundamental lemma of the Calculus of Variations, 
we find that a curve a is a critical point of the action if and only if it satisfies 
5L{a^^[t)) = 0. We have proved the following result. 


Theorem 6.2. An admissible curve aGV{J, E)^^^ is a critical point of the action 
functional S: V{J, E)^^ M. if and only if it satisfies the Euler-Lagrange equations 
5L{af^(t)) = 0 for all tGJ. 
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In order to find the local expression of the Euler-Lagrange equations we just 
calculate the variations of the Lagrangian. To simplify the notation we will use the 
symbol (5y“ for the expression + Cp^y^a^], so that 

A straightforward calculation shows that, for r>2, 




dL 


Sy? 


dt 


r-2 


E(-i) 


d^ f dL 




dP \dy: 




and hence 




_ ^Ol 

Pa^ 


dx^ 

dL 

dx'- 




Jr —1 


r+1 


r—1 


dL 

dP~^ V^2/“ 


_d 

dt 




r=2 j=0 


In the first term, we take into account that <5?/“ = (t“ + and we perform 

a further integration by parts. Denoting by tTo, the coefficient of 6yf, 


(6.7) 

we obtain 


= E(-i) 


r—1 


d^-^ / dL 
dP~^ \dy^ 




_d 

dt 


k r — 2 

r—2 0 


d^ f dL 


dP \dy: 


5Vr-i-l 7 • 


Relabeling the sums we finally get 


{dL,al) = {p, 


dL 


dlTg 

. dt 




( 6 . 8 ) 




E 


k-1 


E(-i)^ 


dp 


dL 


dP- 


dyt 


3 + 1 , 


Sy^r 


Therefore, the Euler-Lagrange equations take the form of the system of ordinary 
differential equations 


(6.9) 


2 ^* = Pgy? 


+ ’"7^^a/32/l 


where the functions tt^ are given by (IgTl) . 
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The coordinate expressions of SL, Tl and 9l can be obtained from the defini¬ 
tion using the local expressions for the vertical endomorphism. However, from the 
integration by parts formula we know that the integrand is { 5L{a?'k) , a) and that 
the boundary terms are (), and hence from equation (16.81) we get 
that the local expression of 5L is 


( 6 . 10 ) 


5L = 


dL 




Taking the induced coordinates (x®, yf,..., yk_i , Mai ■ • ■ i Ma {T^E^)*, the 

local expression of the Legendre transformation El '■ —>■ E^~^)* is of the 

form 


,■ ■ ■, yzfc-i) = ■ ■ ■, yfe_i , Ma,---iMa 

where the are given by the following functions 

k-l 

(6.11) = P^(a^^y?, ■ • •, y2Vi) ^ 


dL 


dy 


1+1, 


Notice that the momenta p)) satisfy the relation 

Po 


ip 1 d Ij , /p 

' =7r^-^TPai 


dy° 


r>l, 


which can serve to define them by recurrence starting from p„ ^ . Also notice 

that p° = tTq,. 

Finally, the coordinate expression of the Cartan form readily follows from that 
of J-L, 


fe-1 


fe-1 


(6.12) y,. = ^p;V“=7r„A“+^p;v; 

with p)) given by (16.1111 . 




r—1 


Noether’s theorem. An immediate consequence of our variational formalism is the 
following version of Noether’s theorem. 


Theorem 6.3. If y is a section of E such that d^kL = dxF for some function 
FgC°°{E^~^), then the function G€C°°{E'^^~^) given by G = F — {9 l ) is 

a first integral. 


Proof. Indeed, we have 

dTG = dTF-dT{9L,ri^^-^) = {dL,ri^)-{dT9L,ri^^) = {SL,r,). 
Therefore, for any solution a of the Euler-Lagrange equations we have 

-f-(G'oa^^“^) = (dTG)oa^^ = { 5 L{a'^^) ,7707) = 0, 
dt 

so that G is constant along any solution of the Euler-Lagrange equations. □ 

In particular, if ry is an infinitesimal symmetry of the Lagrangian, i.e. d^^kL = 0, 
then the function {9l )sC'°“(E^^“^), the momentum in the direction of 77 , 

is a first integral. 
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7. Examples 

We will consider in this section some typical examples of Lie algebroids and we 
will show the form of the Euler-Lagrange equations. 


7.1. The standard case, parameters and quasi-velocities. In the standard 
case E = TM taking a coordinate basis we recover the standard higher- 

order Euler-Lagrange equations [M m uni [3S]- The same equations hold when 
we consider a system defined by a Lagrangian L\{['y]^) = L(A, [ 7 ]^), depending on 
additional parameters AsA. In this case the Lie algebroid \s E = AxTM —j-AxM 
with p(A, v) = (Oa, v) and the bracket is the bracket of vector fields on M depending 
on the variables A as parameters. 

Also, the formalism developed here allows naturally to use a different local basis 
{ei = pI-^} of vector fields in M. In such case the associated coordinates yl are 
called quasi-velocities and the local expressions for the Euler-Lagrange equations 
that we have got is the Euler-Lagrange equations written in quasi-velocities, which 
are sometimes called the higher-order Hammel equations. In that expressions, the 
structure functions are the so called Hammel’s transpositional symbols, defined 
by the equation 


pi- ^ = pic’- 

For more information about quasi-velocities and their use in Mechanics see [29116] 
and for their use in dynamic optimal control see | 2 ]. 


7.2. Systems with holonomic constraints. Let EcTM be an integrable sub¬ 
bundle (i.e. an integrable regular distribution on M). A curve a: R —>■ E is 
admissible if and only if the base curve 7 = roa: R M is contained into an 
integral leaf of E and 0 = 7 . If we denote by E the foliation defined by E, so that 
E = TE, then it follows that = T’^E, that is, it is the set of Ejets of curves 
contained in the leaves of E. 

Given a Lagrangian LgC°°{E’^) an admissible curve a = 7 is a critical point of 
the action if and only if 7 is a critical point of the restriction of the Lagrangian to the 
Etangent T’^{Ej(^tg'j) to the leaf which contains 7 . This follows from standard 

optimization results. Therefore the Euler-Lagrange equations can be obtained as 
the Euler-Lagrange equations for the restriction of the Lagrangian to every leaf, 
usually called the holonomic Euler-Lagrange equations. 

In local coordinates (x®) = adapted to the foliation, so that the leaves 

are given by the equations = constant, we have coordinates 

r = 0,..., fc, in E^ and the Euler-Lagrange equations read 


(7.1) 


k 

E' 

r—O 


-ly 


E dL 


dC dql 


= 0 , 


= k^, 


(r) 


that is, the standard Euler-Lagrange equations on the variables q depending on 
as parameters. 

Alternatively, one can take a local basis {Cq, e^i} of vector fields adapted to the 
distribution E, i.e. E = span({e“}) and then the Euler-Lagrange equations can be 
written as ( 5L{a{t ), Ca ) = 0, where L is any extension of LgC°°{E) to a function 
on TM. In this way we get an expression of the Euler-Lagrange equations with 
holonomic constraints written in terms of quasi-velocities. 
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7.3. Lagrangian systems on Lie algebras. A Lie algebra g can be considered 
as a Lie algebroid over a singleton g {e}. The anchor vanishes from where it 
follows that every curve ^: M —>■ g is admissible and hence we have that g^ is just 
the cartesian product of k copies of g, 

^ (e(0),^(0),^-(0),..., ^(0)) eg X ... X g. 

A section of g —>■ {e} is just an element of g and a local basis {ca} of g provides 
global coordinates (^ 1 ,^ 2 ,...,^*;) in g^. Variational vector fields are of the form 

S^cr = ■ • ■ ’ ^ ^(OJ. 

The Euler-Lagrange equations for a Lagrangian L€C°° {g^) are 
(7.2) TT + adj TT = 0, 

where tt is given by (IgTl) . which in the present case takes the global form 


(7.3) 


TT = 


k 







In this expression ^ stands for the globally defined partial derivative of L with 
respect to given by 


(..., 4 °) , ^) = ^^(4,..., 4-1,4 + sv 4+1,..., 4) 


s=0 


These equations are called the higher-order Euler-Poincare equations [HI [TO] . They 
can be interpreted as the equations for parallel transport of the momentum tt with 
respect to the canonical g-connection on g defined by = [^, ^]. See [33] for the 
details in the first order case. 

The same kind of equations are obtained in the case of a bundle of Lie algebras 
t: K ^ M where the bracket depends on the variables in M, that is, we can have 
different Lie algebra structures on the fibres for different m£M. 


7.4. Systems with advected parameters. We consider a Lie algebra g acting 
on a manifold M by means of a morphism of Lie algebras g X(M); ^ i-+ ^m- 
The trivial bundle r = prj^ : E = M x g ^ M is endowed with a Lie algebroid 
structure with anchor p{m,^) = ^M{rn) and where the bracket is induced naturally 
by the bracket on g (the bracket of constant sections is just the constant section 
corresponding to the bracket on g). 

A curve is admissible if and only if m{t) = It follows 

that they are determined by the curve ^{t) and the initial value rn(0), and hence 
we have the identification E^ = M x g^ given by 

^ (m(0),40),40),..., |3(0)) . 

Variational vector fields are of the form 

^((7 = ^SaTn, Sa^, ■ ■ ■ , — l , 

with Sam{t) = and = a{t) + [40,o’(0]. 
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The Euler-Lagrange equations take the form of the so called Euler-Poincare 
equations with advected parameters 


(7.4) 


TO = ^m(to), 

TT + adj TT = 1 — 1 , 


where tt is given by a expression similar to (1731) but depending also on m, and pm 
is the restriction Pm: 0 —>■ TmM of p to the hbre over m&M, i.e. Pm{0 = ?m(to). 


7.5. Systems with symmetry. We consider a free and proper action of a Lie 
group G on manifold Q, so that p: Q ^ M = Q/G is a principal bundle. The 
vector bundle r = porg: E = TQ/G —>■ M has a natural Lie algebroid structure 
known as the Atiyah algebroid. The anchor is p([w]g) = Tp{v). Sections of E are 
in one to one correspondence to invariant vector fields on Q, and the bracket of 
two invariant vector fields is also an invariant vector field, defining in this way a 
bracket on Sec (A). 

A curve ^(t) in A = TQjG is admissible if and only if there exists a curve q{t) 
on q such that ^{t) = [q{t)]c- It follows that there is a canonical identification 
ijj: T^Q/G {TQ/Gf given by 'ip{[[q]’= ]g) = (where on the left hand 

side the action of G on T^Q is by /c-jet prolongation of the action of G on Q). 

To find some general enough expression for the Euler-Lagrange equations we can 
consider the case of a trivial bundle Q = M x G, so that E = TQ/G = TM xg. 
A curve in E is admissible if and only if v{t) = rhit), where m is the 

base curve of v. Therefore admissible curves are of the form (m(t),^(t)) for m{t) a 
curve in M and ^{t) a curve in g, and we have the identification = T^M x g^, 
given by 

[m,^]'=-i^(H'=;?(0),e(0),...,|S(0)). 

Variational vector fields are of the form 

(C, a) = (c, C ■ ■ •, • ■ ■, > 


with 6cr^{t) = &{t) + [C(t),cr(t)]. 

A Lagrangian L on T^Q/G, generally defined in terms of a G-invariant La- 
grangian on T^Q, produces the following set of Euler-Lagrange equations, which 
are known as the Lagrange-Poincare equations, 

r . dL 

(7.5) 

[71 + ad^ TT = 0 , 

where tt is given by a expression similar to (j7.3|) and ttm is 


ttm 


dL 

5to(i) 





dL 

dmt^k) ■ 


In other words, the equations look like the ordinary Euler-Lagrange equations on 
M together with the Euler-Poincare equations on g. In the case of a non-trivial 
bundle one has to use a principal connection on Q to obtain global expressions for 
the Euler-Lagrange equations (see [mi for the details). 
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8. Reduction and reconstruction 


It was proved in [55] that a morphism of Lie algebroids E ^ E' defines a map 
between the path spaces 3>: V{J,E) —>• V{J,E') by composition $(a) = $oa. If $ 
is fiberwise injective (surjective) then cf> is an immersion (submersion). Moreover, 
variational vector fields are mapped in a simple manner = S$oa(d>otT). 

As a consequence, a morphism induces relations between critical points of functions 
defined on path spaces, in particular between solutions of Euler-Lagrange equations 
for a first-order Lagrangian L in E and a first-order Lagrangian L' = Lo<I> in E', 
producing easy reduction and reconstruction results. 

This correspondence can be easily extended to the higher-order case as follows. 
Consider a Lagrangian LgC°°{E'^) and a Lagrangian L'€C°°{E'^) which are related 
by the map <I>^: E^ that is, L = L' o <I>^. Then the associated action 

functionals S on V{J,E) and S' on V{J,E') are related by d>, that is. S' o ^ = S. 
Indeed, 


S’{^{a)) = S'i^oa) 

rti 

= / (L'o ($ o a)^)(t) dt 

Jto 

rti 

= / {V o o a^){t) dt 

Jto 

ftl 

= / {Loa^){t)dt 

Jto 

= S{a). 


For the boundary conditions, if A'q = ^(Ao) and A'^ = fl>*^ ^(Ai) then $ ap- 

plies 'P{J,E)^^ into V{J, E')j^}, so that the corresponding variational problems 

are related. Indeed, if = Ai then = $^'“^(a^'“^(ti)) = A', 

for i = 0,1. 

Reduction theorems easily follows by considering fiberwise surjective morphisms 
of Lie algebroids. For the sake of clarity, we will omit any reference to the boundary 
conditions. 


Theorem 8.1 (Reduction). Let E ^ E' be a fiberwise surjective morphism of 
Lie algebroids. Consider a Lagrangian L on E^ and a Lagrangian L' on E'^ such 
that L = L' o If a is a solution of the Euler-Lagrange equations for L then 
a' = ^ o a is a solution of Euler-Lagrange equations for L'. 

Proof. Since 5"o$ = S we have that ( d5''(4(a)), Ta ^{ v )) = (d5'(a) , v ) for every 
V € TaV{J, E)^\ If $ is fiberwise surjective, then $ is a submersion, from where 
it follows that $ maps critical points of S into critical points of S', i.e. solutions 
of the Euler-Lagrange equations for L into solutions of Euler-Lagrange equations 
for L'. □ 

We can reduce partially a system and then reduce again the obtained system. 
The final result obviously coincides with the system obtained by the total reduction. 

Theorem 8.2 (Reduction by stages). Let $i: E —>■ E' and ^ 2 ' E' —> E" be 
fiberwise surjective morphisms of Lie algebroids. Let L, L' and L" be Lagrangian 
functions on E^, E'^ and E"^, respectively, such that L'o<I>^ = L and L" = L'. 
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Then the result of reducing first by and later by $2 coincides with the reduction 
6?; <j) = $2 o $1. 

Proof. It is obvious since $ = $2 o $1 is also a fiberwise surjective morphism of 
Lie algebroids. □ 

As an example of the above situation we can consider a system with a group of 
symmetry G. If is closed normal subgroup of G we can reduce first by N and 
later by G/N. Alternatively we can reduce directly by the full group G. The result 
of both procedures is the same. See mm for the first order case. 

For the reconstruction of solutions we have the following result, which is imme¬ 
diate from the variational formalism. 

Theorem 8.3 (Reconstruction). Let E ^ E' be a morphism of Lie algebroids. 
Consider a Lagrangian L on E^ and a Lagrangian L' on E'^ such that L = L'o$^. 
If a is an E-path and a' = ^ o a is a solution of the Euler-Lagrange equations for 
L' then a itself is a solution of the Euler-Lagrange equations for L. 

Proof. Since S' o ^ = S we have that {dS'{^{a)) ,Ta^{v)) = {dS{a) ,v) for 
every v € Ta'P{J,E)^^. If 6 (a) is a solution of Lagrange’s equations for L' then 
dS''( 6 (a)) = 0, from where it follows that d5'(a) =0. □ 

The reconstruction procedure can be understood as follows. Consider a fiberwise 
surjective morphism d>: E ^ E' and the associated reduction map 6 : 'P{J,E) —>• 
P(J, E'). Given an Fl'-path a' G 'P{J, E') solution of the dynamics defined by the 
Lagrangian L', we look for an Fl-path a G P( J, E) solution of the dynamics for the 
Lagrangian L = L' o $, such that a' = 6 (a). For that, it is sufficient to find a 
map T : P( J, E') —>• V{J, E) such that 6 o T = id-pi^j^E')- Indeed, given the Fl'-path 
a' solution for the reduced Lagrangian A', the curve a = T(a') is an A-path and 
satisfy <I>oa = a'. From Reconstruction Theorem 18.31 we deduce that a is a solution 
of the Euler-lagrange equations for the original Lagrangian. 

When specifying the boundary conditions is necessary, the map T must restrict 
toamapTj^J: V{J,E')^} ^V{J,E)jf where AJ, = $'=-i(Ao) and A; = ^'’’-^(Ai). 

Of course, one can define several maps T, and different maps will produce differ¬ 
ent A-paths a for the same A'-path a'. In most cases of interest, $ is fiberwise bijec- 
tive, so that $ is a local diffeomorphism, and the reconstruction process can be done 
with the help of some global diffeomorphism 6 : P( J, E) ^ P x V{J, E'), for some 
manifold P, such that pr 2 06 = 6 . Then fixing pGP we define Tp(a') = 6 “^(p, a'). 
In this way the set of solutions to which a' can be lifted is parametrized by P. 
See bellow for some explicit examples of application of this procedure (and some 
extensions). 

Next, we present some examples where the reduction process indicated above 
can be applied. 

8.1. Euler-Poincare reduction. Consider a Lie group G and its Lie algebra g. 
The map $: TG -G g given by 3>(g,g) = g~^g is a fiberwise bijective morphism 
of Lie algebroids. The fc-jet prolongation : T^G -G g^' of $ is $^([ 5 ]^) = 
T^lg-i[g]^ = [g~^gY'~^, which can be explicitly written as 

<d>\g, 5 ,..., gW) = (e, e, ■ ■ •, with e = g-'g. 
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Assume that L is a left-invariant Lagrangian function on T^G, i.e. 

L{T%-^[g]>^) = L{[h-^g]'^)=L{m 

for every element hGG and every curve g{t) on G. Then we can define a Lagrangian 
L' on by means of 

= Li[g,f) 

where ge is the solution of the differential equation g = Tig^{t) with initial value 
<7e(0) = e. In other words we have L{[g]^) = L'{[g~^g]^~^) for every curve g{t) 
in G, or equivalently L = L'o^^. It follows that if a is a solution of the higher- 
order Euler-Lagrange for L in the Lie group G then a' = $oa is a solution of the 
higher-order Euler-Poincare for L' on g. Moreover, since $ is fiberwise bijective 
every solution can be found in this way, so that the higher-order Euler-Lagrange 
equations on the group reduce to the higher-order Euler-Poincare equations on the 
Lie algebra. 

For the reconstruction process we note that the map : 'P( J, TG) —>■ G x V{J, g) 
given by ^{g,g) = {g{to),T£g-ig) is a global diffeomorphism. Its inverse is the 
map $”^(30,^) = (<?,<?) where g{t) is the integral curve of the left-invariant time- 
dependent vector field X{t,g) = T£g^{t), with initial value g{to) = go. Thus the 
map TgQ(^) = $“^((70,0 provides a reconstruction map. Notice that the curve g{t) 
above is g{t) = goge{t) where ge{t) is the integral curve of X with ge{to) = e. 

8.2. Lie groupoid reduction. For integrable Lie algebroids, the variational prin¬ 
ciple developed here can be obtained via a reduction of an ordinary higher-order 
variational principle with holonomic constraints. In the first-order case this was the 
argument used by Weinstein m to obtain a variational principle on (integrable) 
Lie algebroids. Here we consider the higher-order case. 

Consider a Lie groupoid G over M with source s and target t, and denote by 
E its Lie algebroid, E = A{G). Denote by T^G —?► G the kernel of Ts with the 
structure of Lie algebroid as integrable subbundle of TG. The integral manifolds 
of T^G are the s-fibres of the groupoid, so that T^G is the distribution tangent to 
the foliation <S = UmeMS“^(m), and hence {T^GY = (TS)^ = T^S. 

The map $: T^G —?> E given by left translation to the identities, ‘I*(ug) = 
Tig-i[vg) is a fiberwise bijective morphism of Lie algebroids. As a consequence, if 
L is a Lagrangian function on E^ and L is the associated left invariant Lagrangian 
on (T®G)^, then the solutions of the higher-order Euler-Lagrange equations for L 
(which are the holonomic Euler-Lagrange equations) project by $ to solutions of 
the higher-order Euler-Lagrange equations on the Lie algebroid E. 

For the reconstruction process we consider the manifold G Xm V{J,E) de¬ 
fined by G XmV{J,E) = {{go,a)&G xV{J,E) | t(to) = E(a(to)) }■ The map 
$: 'P(J, T®G) — Gxm'P{J, E) given by $((?,g) = (g(to), Tig-ig) is a global diffeo¬ 
morphism. Its inverse is given by ((70, a) = {g, g) where g{t) is the integral curve 
of the left-invariant time-dependent vector field X(t,g) = Tiga{t) with initial value 
g{to) = e{t{go)) (the identity at the point t{go)). Thus the map TgQ(^) = $“^((70 ,0 
provides a reconstruction map. 

Notice that in this occasion Tg^ is defined as a map Tgg : V{J, E)mo T’( J, T^G)gg 
where T’(J,T^G)go = { {g,g)&'P{J,E) \ g(to)) = go}, mo = t{go) and T’( J, = 

{aGP{J,E) I r(a(to)) = Too }• 
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8.3. Euler-Poincare reduction with advected parameters. Let G be a Lie 

group acting from the right on a manifold M. We consider the Lie algebroid 
E = M X TG ^ M X G where M is a parameter manifold, that is, the anchor 
is p{m,Vg) = (fim^Vg) and the bracket is determined by the standard bracket of 
vector fields on G, i.e. of sections of TG —>■ G, depending on the coordinates in M 
as parameters. Consider also the transformation Lie algebroid E' = M x g ^ M, 
where p{nn,^) = (Cm being the fundamental vector held associated to 

C £ 0). The map $: M x TG M x g given by $(m,Ug) = (mg, g~^Vg) is a 
morphism of Lie algebroids over the action map <p{m, g) = mg, and it is hberwise 
bijective. 

Consider a Lagrangian L on M x T^G depending on the elements of M as 
parameters L{m, [g]^) = Lra{[g]^)- Assume that L is not left invariant but obeys 
to the following transformation rule 

L(m, T^^h[gt) = L{m, [hg]’^) = L{mh, [gf), 

for every element hsG and every curve g{t) in G (equivalently L is invariant by the 
joint left action L{mh~^, [hg]’^) = L{m, [g]’^))- 

We consider the Lagrangian L' on E'^ given by L'(to, = L{m,[ge]^), 

where geit) is the solution of the initial value problem g = g{to) = e. 

Then L' o = L. The variables m which initially are parameters are now dy¬ 
namic variables due to the group action and are called advected parameters. In 
this way, solutions of the higher-order Euler-Lagrange equations for L (standard 
Euler-Lagrange equations with parameters) are mapped by $ to solutions of the 
higher-order Euler-Lagrange equations for L', i.e. the higher-order Euler-Poincare 
equations with advected parameters. 

For the reconstruction process we consider the global diffeomorphism $: V{J,Mx 
TG ) -)> G X V{J,M X 0) given by $(m, (5,5)) = (^(to), { mg , Tlg ~ ig )) . Its inverse 
is (w,C)) = (5i5)) where g[t) is the integral curve of the left- 

invariant time-dependent vector field X{t,g) = Tig^{t) with g(to) = go- Thus the 
map Tgg{m,^) = ^~^[go, (m,C)) provides a reconstruction map. 

8.4. Lagrange-Poincare reduction. We consider a Lie group G acting free and 
properly on a manifold Q, so that the quotient map p: Q ^ M = Q/G has the 
structure of a principal bundle. We consider the standard Lie algebroid structure 
on E = TQ and the Atiyah algebroid E' = TQ/G —>• M. The quotient map 
$: TQ —)• TQ/G, <I>(u) = is a Lie algebroid morphism and it is hberwise 
bijective. Every G-invariant Lagrangian on T^Q dehnes uniquely a Lagrangian L' 
on E'^ such that L' o = L, explicitly given by 

L'mG]^-^) = L'{[[q]^]G) = L{[q'^]). 

Therefore every solution of the G-invariant Lagrangian on T^Q projects to a solu¬ 
tion of the reduced Lagrangian on {TQ/G)^ = T^Q/G, and every solution on the 
reduced space can be obtained in this way. Thus, the higher-order Euler-Lagrange 
equations on the principal bundle reduce to the higher-order Lagrange-Poincare 
equations on the Atiyah algebroid m- 

For the reconstruction process we consider the manifold Q Xm V{J,TQ/G) 
dehned hy Q Xm 'P{J,TQ/G) = {{qo,a)GQ x'P{J,TQ/G) | p(to) = r(a(to)) }■ 
The map $: V{J,TQ) ->■ Q Xm v\j,TQ/G) given by $(g,d) = (9(^0), [die) is 
a global diffeomorphism. Its inverse is given by $“^(go,o) = {g^g) where q{t) 
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is the curve determined as follows. The curve a{t) is admissible, so that it is of 
the form a{t) = [q{t), q(t)]c', if go€G is the unique element in the group such 
that q{to) = goqo, then we define q{t) = g^^q{t). Therefore the map Tqg{a) = 
^~^{qo,a) provides a reconstruction map, which in this case is defined as a map 
Tgg : V{J,TQ/G)mo 'P{J,TQ)qg where mo = p{qo). 

9. Conclusions and outlook 

In this paper we have obtained the Euler-Lagrange equations for a higher-order 
Lagrangian. The emphasis has been on the variational description using the tools 
of variational calculus on Lie algebroids developed in [5S] . 

However, there are other many aspects of the theory that we have left out 
and will be studied in future. The variational structure strongly suggests that 
a symplectic formalism is possible, similar to Klein’s formalism |25j of standard 
Lagrangian Mechanics and its generalization to Lie algebroids [2^, and the corre¬ 
sponding Hamiltonian version [301 ng. It is interesting to study such formalisms 
and the transformation properties of the symplectic form. 

On the other hand, a geometric version of Pontryagin maximum principle which 
allows reduction in the setting of Lie algebroids was studied in Eiiisa (a direct 
proof was given in [23], based on needle variations [38] and the results in [3^). It 
is interesting to study the relation between that results with those on this paper 
and the proposed Hamiltonian versions already mentioned, on one hand, and with 
other ‘variational principles’ that appeared recently in the literature (Hamilton- 
Pontryagin [55], Clebsh-Pontryagin m, etc.), on the other. 

Concerning the optimality properties of the solutions, in this paper we have ob¬ 
tained only first order conditions, that is conditions for critical points of the action. 
To characterize those which are local maxima/minima we will study further opti¬ 
mality conditions. Since our results are based on a bona fide variational principle 
it is expected that the Hessian of the action contains such information. Moreover, 
when two Lagrangians are related by a morphism the Hessians at corresponding 
critical points should be also related. 

Problems with time dependent Lagrangians can be treated by embedding E into 
TM. X E as it is mentioned in [36] . However the formalism introduced in [SS] [22] for 
first-order systems can be easily generalized to the present case. 

All this problems and the possible generalizations to field theory, following the 
ideas in [33], will be studied elsewhere. 

References 

[1] Abrunheiro L and Camarinha L 

optimal control of affine connection control systems from the point of view of lie algebroids 
Int. J. Geom. Methods Mod. Phys. (2014) DOI; 10.1142/S0219887814500388 

[2] Abrunheiro L, Camarinha M, Carinena JF, Clemente-Gallardo J, Martinez E and 
Santos P 

Some applications of quasi-velocities in optimal control 

International Journal of Geometric Methods in Modern Physics,8 04, 2011, 835-851 

[3] Camarinha M, Silva-Leite F and Crouch P 
Splines of class on non-Euclidean spaces 
IMA J. Math. Control Info. 12 (4), 1995, 399-410 

[4] Carinena JF, Lopez C and Martinez E Sections along a map applied to higher-order 
Mechanics. Noether’s Theorem 

Acta Applicandae Mathematicae 25 (1991) 127-151 


HIGHER-ORDER VARIATIONAL CALCULUS ON LIE ALGEBROIDS 


27 


[5] Carinena JF and Martinez E 

Lie Algebroid generalization of Geometrig Mechanics 

In Lie Algebroids and related topics in differential geometry, Banach Center Publications 54, 
2001, p. 201 

[6] Carinena JF, Nunes da Costa JM and Santos P 
Quasi-coordinates from the point of view of Lie algebroid structures 
Journal of Physics A: Mathematical and Theoretical, 40, 2007, 10031-10048 

[7] Cendra H, Marsden JE and Ratiu TS 
Lagrangian reduction by stages 

Mem. Amer. Math. Soc. 152 (2001), no. 722, x+108 pp 

[8] Colombo L 

Lagrange-Poincare reduction for optimal control of underactuated mechanical systems 
Preprint available at http://arxiv.org/pdf/1306.6005 

[9] Colombo L and Martin de Diego D 

A variational and geometric approach for the second-order Euler-Poincare equations 
Notes of the talk delivered at XIII Encuentro de Invierno, Zaragoza (2011) 
available online at http://andres.unizar.es/~ei/2011/Contribuciones/LeoColombo.pdf. 

[10] Colombo L and Martin de Diego D 

On the geometry of higher-order variational problems on Lie groups 
Preprint available at http://arxiv.org/pdf/1104.3221vl. 

[11] Colombo L, Martin De Diego D and Zuccalli M 

Optimal control of underactuated mechanical systems: a geometric approach 
J. Math. Phys. 51 (8), 2010, 083519 

[12] Cortes J, de Leon M, Marrero JC and Martinez E 
Nonholonomic Lagrangian systems on Lie algebroids 

Discrete and Continuous Dynamical Systems A 24 (2), 2009, 213—271 

[13] Crainic M and Fernandes RL 
Integrability of Lie brackets 

Ann. of Math. (2) 157 (2003), no. 2, 575-620 

[14] Crampin M, Sarlet W and Cantrijn F 

Higher order differential equations and higher order Lagrangian mechanics 
Math. Proc. Camb. Phil. Soc. 99 (1986) 565. 

[15] de Leon M, Marrero JC and Martinez E 
Lagrangian submanifolds and dynamics on Lie algebroids 
J. Phys. A: Math. Gen. 38 (2005), R241-R308 

[16] De Leon M and Rodrigues P 

Generalized Classical Mechanics and Field Theory 
North-Holland Math. Stu. 112, Amsterdam, 1985. 

[17] Gay-Balmaz F, Holm DD, and Ratiu TS 

Higher order Lagrange-Poincare and Hamilton-Poincare reductions 
J. Braz. Math. Soc. 42(4), (2011), 579-606 

[18] Gay-Balmaz F, Holm DD, Meier DM, Ratiu TS and Vialard FX 
Invariant Higher-Order Variational Problems 

Comm. Math. Phys. 309 (2012), 413-458. 

[19] Gay-Balmaz F and Ratiu TS 

Clebsch optimal control formulation in mechanics 
J. Geometric Mechanics. 3 (1), (2011), 41-79 

[20] Grabowska K and Grabowski j 

Variational calculus with constraints on general algebroids 
J. Phys. A: Math. Theor. 41 (2008), 175204. 

[21] Grabowska K, Grabowski J and Urbanski P 
Geometric mechanics on algebroids 

Int. J. Geom. Meth. Mod. Phys. 3 (2006), 559—575. 

[22] Grabowska K, Grabowski J and Urbanski P Lie brackets on affine bundles Ann. Global 
Anal. Geom. 24 (2003), 101-130 

[23] Grabowski J and Jozwikowski M 

Pontryagin maximum principle on almost Lie algebroids 
SIAM J. Control Optim. 49 no. 3, 2011, 1306-1357 


28 


E. MARTINEZ 


[24] JOZWIKOWSKI M AND ROTKIEWICZ M 

Prototypes of higher algebroids with applications to variational calculus 
Preprint available at http://arxiv.org/abs/1306.3379 

[25] Klein J 

Espaces variationnels et mecanique 
Ann. Inst. Fourier 12 (1962) 1-124 

[26] Kyriakopoulos KJ and Saridis GN 
Minimum jerk path generation 

In: Proceedings of the 1988 IEEE International Conference on Robotics and Automation, 
Philadelphia, 1988, pp. 364-369 

[27] Machado L, Silva-Leite F, and Krakowski K 

Higher-order smoothing splines versus least squares problems on Riemannian manifolds 
J. Dyn. and Control Syst. 16, 2010, 121-148 

[28] Mackenzie KCH 

General Theory of Lie Groupoids and Lie Algebroids 
Cambridge University Press, 2005 

[29] Martinez E 

Lagrangian Mechanics on Lie algebroids 
Acta Appl. Math., 67 (2001), 295-320 

[30] Martinez E 

Geometric formulation of Mechanics on Lie algebroids 

In Proceedings of the VIII Fall Workshop on Geometry and Physics, Medina del Campo, 
1999, Publicaciones de la RSME, 2 (2001), 209-222 

[31] Martinez E 

Reduction in optimal control theory 
Rep. Math. Phys. 53 (2004) 79-90 

[32] Martinez E 

Lie Algebroids in Classical Mechanics and Optimal Control 
SIGMA 3 (2007), 50, 17 pages 

[33] Martinez E 

Classical Field Theory on Lie algebroids: Variational aspects 
J. Phys. A: Mat. Gen. 38 (2005) 7145-7160 

[34] Martinez E 

The momentum equation 

In Groups, Geometry and Physics. Monografias de la Real Academia de Ciencias de Zaragoza 
29 (2006) 

[35] Martinez E, Mestdag T and Sarlet W Lie algebroid structures and Lagrangian systems 
on afBne bundles J. Geom. Phys. 44 (2002), no. 1, 70-95 

[36] Martinez E 

Variational calculus on Lie algebroids 

ESAIM: Control, Optimisation and Calculus of Variations 14 02, 2007, 356-380 

[37] Noakes L, Heinzinger G and Paden B 
Cubic splines on curved spaces 

IMA Journal of Mathematical Control &; Information 6 1989, 465-473 

[38] PONTRYAGIN LS, BOLTYANSKII VG, GAMKRELIDZE RV Y MiSHCHENKO EF 
The mathematical theory of optimal processes 

Interscience Publishers John Wiley &: Sons, Inc. New York-London, 1962 

[39] Prieto-Martinez PD and Roman-Roy N 

Lagrangian-Hamiltonian unified formalism for autonomous higher-order dynamical systems 
J. Phys. A: Math. Theor. 44 (2011), 385203. 

[40] Tulczyjew W 

The Lagrange differential 

Bull. Acad. Polon. Sci. 24 (1976), 1089-1097. 

[41] Weinstein A 

Lagrangian Mechanics and groupoids 
Fields Inst. Comm. 7 (1996), 207-231 


HIGHER-ORDER VARIATIONAL CALCULUS ON LIE ALGEBROIDS 


29 


[42] Yoshimura H and Marsden JE 

Reduction of Dirac structures and the Hamilton-Pontryagin principle 
Rep. Math. Phys. 60, 2007, 381-426 

Eduardo Martinez: IUMA-Departamento de Matematica Aplicada, Universidad de 
Zaragoza, Pedro Cerbuna 12, 50009 Zaragoza, Spain 
E-mail address: emfSposta.uiiizar.es 


